Rapid Merger of Binary Primordial Black Holes 



o 
o 

(N 

in 



Kimitake Hayasaki^^^jB Kcitaro Takahashi^-'^, Yuuiti Scndouda^, and Shigchiro Nagataki^ 
^ Yukawa Institute for Theoretical Physics, Kyoto University, 
Kita-shirakawa Oiwake-cho, Sakyo-ku, Kyoto 606-8502, Japan 
^Division of Physics, Graduate School of Science, 
Hokkaido University, Kitaku, Sapporo 060-0810, Japan and 
'^Department of Physics, Nagoya University, Chikusa-ku, Nagoya 4.6^-8602, Japan 

(Dated: September 9, 2009) 

We propose a new scenario for the evolution of a binary of primordial black holes (PBHs). We 
consider the dynamical friction by ambient dark matter and gravitational interaction between a 
binary and a circumbinary disk, assuming PBHs do not constitute the bulk of dark matter. After 
the turnaround, a PBH binary loses the energy and angular momentum by the two processes, 
which are very effective for a typical configuration. Finally the binary coalesces due to the emission 
of gravitational waves in a time scale much shorter than the age of the universe. We estimate 
the density parameter of the resultant gravitational wave background. Astrophysical implication 
concerning supermassive black holes is also discussed. 

PACS numbers: 95.85.Sz,98.62. Js,98.80.Cq 
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Introduction. Primordial black holes (PBHs) are hy- 
pothetical objects formed due to gravitational collapse of 
density perturbations at the early phase of the universe 
A PBH mass larger than ~ 10"1^Mq, where Mq is 
the solar mass, is enough to prohibit evaporation via the 
Hawking radiation and PBHs with such a mass would 
remain as massive compact halo objects (MACHOs) un- 
til now. The possibility that PBHs consitute a certain 
fraction of dark matter has been investigated by many 
authors (see 0] and references therein). 

In Q, Nakamura et al. considered the formation and 
evolution of PBH binaries. If PBHs were formed ran- 
domly rather than uniformly in space, some pairs of 
PBHs would have sufficiently small separations to over- 
come the cosmological expansion and then form a bound 
system, binary. They estimated gravitational waves from 
the coalescence of the binaries following the binary evo- 
lution taking gravitational wave emission into account. 
Because gravitational radiation is not effective to extract 
angular momentum from a binary unless the separation 
is very small, only a small fraction of binaries coalesces 
during the age of the universe. Nontheless, cosmological 
gravitational wave background amounts to Slgw ~ 10" 
at the LISA band (i^gw ~ 10""* Hz), where flgw and i/gw 
are the density parameter and frequency of gravitational 
waves, respectively (see also 0)- 

In this Letter, we consider two processes which affect 
the evolution of a PBH binary. One is the dynamical 
friction by ambient dark matter and another is the inter- 
action between a binary and a gaseous disk surrounding 
it (i.e. circumbinary disk). Both processes are inevitable 
and help the binary shrink. Our scenario predicts that 
PBH binaries generically coalesce during the age of the 
universe and gravitational waves are emitted more effi- 
ciently compared to the result in Throughout this 
paper, we adopt the cosmological parameters obtained 
by 5-yr WMAP 



Binary PBH. First we give a summary on the evolution 
scenario of PBH binary presented by [1, 0] , extending to 
a case where PBHs do not consitute the bulk of dark 
matter. For simplicity, we assume a primordial spectrum 
of density fluctuations which has a sharp peak with a 
sufficient amplitude to form PBHs so that they have a 
monochromatic mass function (see, e.g., [lo| for such an 
inflationary model). The radiation and matter inside a 
horizon gravitationally collapse into a black hole if the 
over density is large enough there. The black hole mass, 
A'/bh, is comparable to the horizon mass and can be writ- 
ten as Mbh ~ (Tbhf/l GeV)-2MQ, where IbM is the 
temperature of the universe at the formation [ll|. We 
assume that PBHs constitute a fraction, /, of the dark 
matter, that is, fibh = /f^DM where fibh and JIdm are 
the current density parameter of PBHs and dark mat- 
ter, respectively. The value of / is constrained for a wide 
range of Mbh and can be as large as 0.1 for Mbh O.IMq 
while much severer constraints, / < 10~^, are obtained 
for Mbh > IO^Mq d. 

Noting that the density of PBHs can be written as 
Pbh(T) = /pdm(T) = /TeqT3 where Teq ~ 1 eV is the 
temperature at the matter-radiation equality, the average 
separation of PBHs at the formation, f, is expressed as 
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Actually, PBHs are randomly distributed in space at the 
formation and we consider a pair of nearest PBHs with 
the initial separation parametrized as af. Statistically, 
most of nearest PBH pairs will have a < 1 while only a 
few would have a ^ 1. In fact, the probability density 



for a random 



distribution function is dP/da 
distribution as given by 

After the black hole formation, the separation of the 
pair evolves in proportion to the scale factor due to the 
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cosmic expansion. However, if the local energy density of 
the pair, ppair = <^^^Phh, becomes larger than the radia- 
tion energy density, p^, the pair decouples from the cos- 
mic expansion and forms a gravitationally bound system. 
We see that the turnaround occurs at the temperature. 
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Here it should be noted that for a > f^^^ we have always 
Ppair < Pdm SO that tlic turnaround never occurs. Thus 
wc obtain a condition for a that a pair forms a bound sys- 
tem, ao < a < where ao = f^/^M~^''^Tll^ M}''^ 
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2 X 10-3/V3(j\^j^j^/ijVfQ)i/6 ^Yie lower bound comes 
from the requirement that af should be greater than the 
Schwarzschild radius of a PBH. The separation at the 
turnaround is. 
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After the turnaround, if the two black holes had no rel- 
ative velocity, they would coalesce to form a single black 
hole. However, as discussed in [^Q, the tidal force from 
neighboring black holes would give the pair some angu- 
lar momentum and prevent the coalescence. From the 
above arguments, the typical separation between bina- 
ries is given by a /^/'^ while another PBH which gives 
angular momentum to the binary would be typically sep- 
arated from the binary by ~ f. The semi-major axis of 
the binary is roughly given by a ~ rta. On the other 
hand, the acceleration due to the tidal force and free fall 
time are estimated as, 
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respectively, where R = (Tbhf/Ttajr is the typical sep- 
aration between the binary and the third PBH at the 
turnaround of the pair. Thus the semi-minor axis is es- 
timated as & ~ fltidai^ff ^ '''ial^^- We Can see that the 
typical eccentricity of the binary is e = -y/l — b'^ / a? ~ 

Once a binary is formed, it emits gravitational waves 
and the separation shrinks gradually. The coalescence 
timescale of a binary PBH is given by [l^ . 
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Thus we see that the emission of gravitational waves is 
not effective for most of the binaries {a < f'^^^) to merge 



during the age of the universe, unless / is very small 
/ < 10-3(Mbh/lMQ)5/2i^ However, as we see below, 
considering two other processes which extract the energy 
and angular momentum from the binary, its evolution 
drastically changes. Hereafter we consider a typical pair 
with a < Note that if / is very small, the number 

of PBH binaries is very small so that the total amount 
of gravitational waves is also small even if most of PBH 
binaries merge during the age of the universe. 

Interaction with dark matter. In the case of / ^ 1 we 
consider here, dark matter other than PBHs is abundant 
and will accrete onto PBHs forming halos around them. 
The separation of a binary in a halo decays due to the 
dynamical friction from dark matter. Below we argue 
how effectively this happens. 

At the turnaround of a pair of PBHs, each PBH would 
have a halo with the density ~ pumizta) and mass, 
Mhaio ~ PDM^ta ^ (a^//)A-^bh, which is comparable 
to the PBH mass. It would be reasonable to assume 
that, at the turnaround, two halos merge into a single 
halo which surrounds the binary. Further, ambient dark 
matter continues to accrete onto the binary. In matter 
dominant era, the halo mass and turnaround radius of 
dark matter evolve as Afhaio ^ A/bh(l + Zcq)/i^ + z) and 

rta.dm ~ lO-2{Mbh/lM0)l/3((l + ^)/(l + z,q))-4/3 

respectively [13|, where z^q ~ 3000 is the redshift at the 
matter-radiation equality. 

A binary in a halo shrinks by giving angular momen- 
tum and energy to dark matter particles. However, the 
binary-dark matter interaction can happen only if dark 
matter particles pass enough close to the binary, in other 
words, if they enter the loss cone Also, how much 
the binary shrinks depends on the dark matter mass 
which enter the loss cone. Here we simply assume that a 
dark matter particle enters the loss cone if its periccntcr 
is smaller than the semi-major axis of the binary. Let 
us evaluate the orbital elements of a typical dark mat- 
ter particle. As in the case of a PBH pair, dark matter 
receives tidal force from nearby PBHs. The semi-major 
axis and semi-minor axis are again estimated as ~ rta, dm 



and 
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ta,dm -, „ — 2 f2 



i?6 



1 + z \-^"/^ /ATbh 



1 + z. 



cq 



1/3 



pc. 



From d < rta, wc obtain the following condition: 
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This means that dark matter which turnarounds before 
this redshift enters the loss cone. The total mass which 
enters the loss cone is then bounded by, 

„3 \ 2/5 

(8) 
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Next we consider the time scale of the dynamical fric- 
tion. It depends on the dark matter density of the halo 
and we assume it to be equal to the average density of 
the universe. In fact, the dark matter would be denser in 
a halo and this assumption overestimates the timescale, 
which leads to a conservative result. The time scale is 
given by 
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where erf is the error function, In A « 10 is the Coulomb 
logarithm. Here X = Ubh/(V2cr) where Wbh and a are 
the black hole velocity and velocity dispersion of dark 
matter, respectively. We conservatively assume that the 
dynamical friction is effective only when this time scale 
is smaller than the Hubble time at that time. In matter 
dominated era, the ratio is, 
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where F '--^ 20 is a numerical factor and we set and 
a to the Kcplcrian velocity y^^Mbh/r^. Thus the dy- 
namical friction is effective for 1 + z > r~^/'^(l -I- Zoq). 
The total mass which has entered the loss cone by this 
time is estimated as Mic ^ F^/'^Afbh- Together with the 
argument around Eq. ([8]), our estimation of the total 
mass which contributes to the dynamical friction is Mic ~ 
min[/-3/5Af|^j^,r2/3Afbh], which is T^/^A/bh 7Afbh for 
/ ^ 0.04. 

Once we obtain the loss cone mass, we can evaluate 
how much the binary shrinks, using the fact that the 
mean energy which dark matter of mass dA/ic extracts 
from a binary is ^ GAI\j\idM\c/ a: 
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Thus, the resultant separation after the dynamical fric- 
tion ceases, a^f, is given by. 
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Adopting a typical value given above, M\c ^ 7Ajfbh, we 
see that the binary shrinks by about three orders by z ~ 
400. Here it should be noted that, in the above analysis, 
we used several simplifying assumptions. In particular, 
there would be a large uncertainty in the estimation of 
the shrinking factor. Nevertheless, we believe that the 
dynamical friction would have an important effect on the 
evolution of PBH binaries unless / is close to unity. More 
robust evaluation of the shrinking factor would need N- 
body simulations, which we will present elsewhere. 

Interaction with gas. Baryon g well as the dark 

matter, also accretes onto the PBHs. Even if the gas has 



negligible angular momentum with respect to the cen- 
ter of mass of the binary, the gas can not fall directly 
into the black holes but instead would form a rotating 
disk around the binary (i.e. a circumbinary disk). The 
circumbinary disk can play a role of an efficient mech- 
anism to extract the angular momentum from a binary 
16l . 17 , [H, [l^ . This is because the circumbinary disk 
and binary exchange their masses and ang ular momenta 
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through the tidal-resonant interaction 

Neglecting the angular momentum of the gas, we con- 
sider the Bondi accretion and its accretion radius is given 
by 
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where Coo~ 10(l + z)/(l + Zcq) km/s is the average sound 
velocity [23[ . When the Bondi radius is smaller than the 
separation of the binary, the gas accretes onto each PBH 
separately and a circumbinary disk would not form. Be- 
cause the Bondi radius increases in time while the binary 
separation would shrink due to the dynamical friction, 
the Bondi radius exceeds the separation sooner or later 
and a circumbinary disk would form. If / is relatively 
small / ^ 0.1 and the dynamical friction is as effective 
as we argued in the previous section, a circumbinary disk 
would begin to form at z < 1000. 

According to [1, [2^, the gas accretion rate onto a 
black hole with 1A/q has a peak around the matter- 
radiation equality time and the peak accretion rate is 
several percent of the Eddington rate: A/Edd ~ 2 x 
lO"9(Afbh/lM0)Afo/yr. For A/bh < IMq, the mass of 
the circumbinary disk Afcbd ~ -^Edd^^"^ can be well fit- 
ted by A/cbdA/bh ^ 10-5(Afbh/lMo). Considering only 
the interaction between the binary and the circumbinary 
disk, that is, neglecting the effects of the respective disks, 
the time scale of the orbital decay can be estimated as 
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where tvis ~ 7 x 10'^(a/lpc)^/^yr is the viscous time scale 
of the circumbinary disk with an assumption that the 
temperature of the inner edge of the circumbinary disk is 
(1 + z)Tq and the Shakura-Sunyaev viscosity parameter is 
0.1 [2J|. Further we used a typical parameter set, a^f ^ 
10^"* pc and 1 + z ^ 400 taking the dynamical friction 
and time evolution of the Bondi radius into account. We 
will use this parameter set below as well. 

As the separation decreases, the timescale ([T4)) also 
decreases gradually while the timescale of gravitational 
waves, Eq. ([5|), decreases much more rapidly. When the 
two time scales become comparable, the most effective 
process of extracting angular momentum from the binary 
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switches to the gravitational wave emission. The critical 
separation, Oc, is, 
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and the critical redshift, z^, which is determined by 
H^^ = tcbd is given by. 
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Here we assumed the binary orbit is circular after the 
dynamical friction. Substituting Eq. (fT5|) into Eq. ((5]), 
the time scale of gravitational radiation reduces to ig^ ^ 
10'*(Mbh/lMo)-i^/2iyj.^ Thus, for the typical parameter 
set, PBH binaries of Mbh ^ O.OIMq generically coalesce 
by around Zc > 5. This is to be contrasted with Q where 
only a small fraction of binary PBHs coalesces during the 
age of the universe. 

Summary and Discussions. We have studied the evo- 
lution of PBH binaries by considering the gravitational 
interactions with dark matter and baryon gas. Let us 
summarize our scenario in turn. After the turnaround 
of a pair of PBHs {z > ^oq), dark matter accretes onto 
the pair to form a halo. In the early phase (z < ^cq), 
the dark matter density is sufficiently high that the dy- 
namical friction is effective. Although a precise estima- 
tion of the shrinking factor of the binary separation is 
hard analytically, our simple model implies that a binary 
can shrink by a couple of orders. Some time after the 
recombination (z < 1000), a circumbinary disk would 
form and extract angular momentum from the binary. 
At z ~ O(IOO), the most effective process extracting an- 
gular momentum from the binary switches to the gravita- 
tional wave emission. Finally the binary would coalesce 
at around z ~ 100. 

If our scenario works, most of PBH binaries coalesce 
around z = Zc emitting gravitational waves. Thus we can 
expect a substantial amount of the gravitational wave 
background. Below we present an order-of-magnitude 
estimation simply assuming that all binaries coalesce at 
z — Zc- Given the probability distribution function for 
the initial separation, dP/da ^ 3a^e~" , we see that the 
fraction of PBHs which form a binary is P{a < f^^^) ^ f 
assuming / <JC 1. Then, at z — Zc, the number of bina- 
ries in a horizon volume is given by ~ pH~^ pj^M./M\,\i. 
Noting that the energy of gravitational waves emitted by 
a binary is ~ (GAfbhJ^)^''^Afbh, where v = y^GMbh/a^ is 
the frequency of the gravitational wave determined by the 
binary separation, the energy density of the gravitational 
waves at z = Zc is e^^{v,Zc) ~ {GM\,)^vY/^ p pj^m{zc)- 
Finally, the current density parameter of the gravita- 



tional waves is given by, 
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where z/c,o - 5 x lO-\M^^/MQf/^\a/a,)-^/^/{l + 
Zc)Hz represents the lowest frequency corresponding to 
the separation Oc- The frequency domain is given as 
Vcfl < ly < t'mso,o where 105(Mbh/Mo)-V(l + 

Zc)Hz, which is the frequency corresponding to the 
marginally stable orbit of a PBH binary. The gravita- 
tional wave background with / < 0.03 is estimated as 
^^gw,o ~ 10"^*^, which is comparable with that estima- 
tion by corresponding to / = 1, and is enough large 
to be detected with LISA, DECIGO and LIGO. Thus the 
observation of the gravitational wave background would 
give a strong constraint on the dark matter fraction of 
PBHs with Mbh < ^Mq. 

Finally let us mention another astrophysical implica- 
tion of our scenario. If the spectrum of density pertur- 
bation is broad rather than monochromatic as assumed 
here, clusters of PBHs can be formed [2^. If this is the 
case, PBH mergers may occur frequently in clusters to 
form much larger black holes than the original PBHs. 
This may help explain the presence of a supermassive 
black hole recently discovered at z ~ 6 2^, 27 1, providing 
supermassive black holes directly by successive merger of 
PBHs or smaller "seed" black holes. 
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